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We study the degeneracies between neutrino mass and dark energy as they manifest themselves in 
cosmological observations. In contradiction to a popular formula in the literature, the suppression 
of the matter power spectrum caused by massive neutrinos is not just a function of the ratio of 
neutrino to total mass densities — fl^/flm, but also each of the densities independently. We 
also present a fitting formula for the logarithmic growth factor of perturbations in a flat universe, 
/(z, fc; w, Qde) ^ [1 — A{k)Q,DEfv + B{k)fl — C {k) f^]Q'^{z) , where a depends on the dark 
energy equation of state parameter w. We then discuss cosmological probes where the / factor 
directly appears: peculiar velocities, redshift distortion and the Intergrated Sachs- Wolfe effect. We 
also modify the approximation of Eisenstein & Hu (1999) for the power spectrum of fluctuations in 
the presence of massive neutrinos and provide a revised cod^. 



I. INTRODUCTION 

The latest results from the WMAP satellite [l| con- 
firm the success of the ACDM model, where ~ 75 % of 
the mass-energy density is in the form of dark energy, 
and matter, most of it in the form of Cold Dark Matter 
(CDM) making up the remaining 25 %. Neutrinos with 
masses on the eV scale or below will be a hot component 
of the dark matter and will free-stream out of overdensi- 
ties and thus wipe out small-scale structures. This fact 
makes it possible to use observations of the clustering of 
matter in the universe to put upper bounds on the neu- 
trino masses. A thorough review of the subject is found 
in 0]. With the improved quality of cosmological data 
seen in recent years, the upper limits have improved, and 
some quite impressive claims have been made in the re- 
cent literature, e.g. 

Present cosmological neutrino mass limits make use of 
the suppression effect of the neutrino free-streaming at a 
fixed, given redshift. As our ability to map out the mass 
distribution at different epochs of the cosmic history im- 
proves, by doing, e.g., weak lensing tomography, we will 
gain sensitivity by in addition using the effect of massive 
neutrinos on the growth rate of density fluctuations. One 
key issue which then arises is possible degeneracies be- 
tween neutrino masses and cosmological parameters. In 
this paper we focus on the degeneracy between the dark 
energy equation of state and the neutrino masses. We 
show that the combined effect of neutrinos and dark en- 
ergy can be parametrized in a simple manner, and that 
the degeneracy can be broken by mapping out the large- 
scale structure over a reasonably wide range of redshifts. 
We also explore the effect of massive neutrinos on the 
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TABLE I: The ability of observational tests based on geome- 
try, matter power spectrum and linear growth to probe neu- 
trino mass and dark energy. In the text we discuss the sensi- 
tivity of P(k, z) to neutrino mass. 



matter power spectrum P(k), and we modify the formula 
given by Eisenstein & Hu [4] , such that it would be valid 
for realistic properties of massive neutinos. The depen- 
dencies of three important ingredients of the universe on 
the equation of state parameter w and the neutrino den- 
sity Vli, are shown in Table IH 

The outline of the paper is as follows. In Section II we 
contrast common approximations to the power spectrum 
of fluctuations with the exact results from CAMB, and 
we provide a new modified approximation. In Section 
III we provide a new fitting formula for the linear the- 
ory growth of perturbations in the presence of massive 
neutrinos. In section IV we discuss the parameter de- 
generacy between dark energy parameters and neutrino 
masses, and we illustrate how it manifests itself in pecu- 
liar velocities and the Integrated Sachs Wolfe effect. Our 
conclusions are summarized in Section V. 



II. MASSIVE NEUTRINOS AND THE MATTER 
POWER SPECTRUM 

Most cosmological neutrino mass limits make use 
of the matter power spectrum P{k) in some guise, 
although it is also possible to obtain a limit from cosmic 
microwave background data alone, see [1, d, 0|- 
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A useful way to consider the effect of neutrino mass on 
the power spectrum is to consider the following quantity 
at a given redshift: 



AP(fc) Pik; U) - P{k; U = 0) 



P{k) 



P{k; U = 0) 



(1) 



A common heuristic explanation for the role of the mat- 
ter power spectrum in deriving neutrino mass limits is 
the approximate expression « — 8/i/, describing 

the suppression of small-scale power caused by neutrino 
free-streaming, a result valid for 1, and first given in 
Q. See [1, Q for derivations. In Figure ([T]), we compare 
this approximation to the matter power spectrum using 
C AMB 1 [10] , but also to the Eisenstein & Hu approxima- 
tion to our modification to this E&H approximation 
and to the numerical solution given by equation for 
fi, ~ — 0.04 and = 0.16. An approximation for 
P(fc, z) is given in [2] 



P(fc,^) = < 



(l + z)5(0), 
(1 + 2)5(0) 



P(fc,0) ioraH<k<kn 



P{k,Q) for 



k-» k„ 



(2) 

where p ~ 1 — 3/5/^ based on equation ([7]) below. We 
give a better approximation to this factor in Section FlIIBI 
The non-relativistic scale, fc„r, is given by 



0.01817^/2 (-^)i/2 h Mpc- 



(3) 



For another approximation valid for Mixed Dark Matter 
see [13. 

As a rule of thumb, the present-epoch matter power 
spectrum is considered to be in the linear regime for 
comoving wavenumbers k < 0.10-0.15 hMpc^^, and we 
see from Figure (P) that AP/P obtained from CAME 
tends to a constant only for k > 1.5 Mpc^^, which is 
well into the non-linear regime of structure formation. 
This fact is also evident from Figures (12) and (13) 
in 0. Thus, ^ « Sfi, should only be used as a 
heuristic guide to the effect of massive neutrinos on 
the power spectrum since it's not valid at the length 
scales where the matter power spectrum can be said to 
be in the linear regime. Furthermore, we note that, as 
expected, it works well only for very small fi, whereas for 
large neutrino masses this approximation breaks down. 
Moreover, Figure ^ also shows that for small neutrino 
masses the E&H approximation breaks down. However, 
by modifying the master transfer function used by E&H, 
we managed to minimise the error between CAMB and 
the E&H approximation for = 0.04, from 20% to only 
3% for 0.02 < A: < 0.15 h Mpc"\ This modification 
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FIG. 1: The quantity AP/P a.t z — defined in equation 
(HI, derived using CAMB for fixed w = -1, = 0.25, fib = 
0.04, h = 0.7 (solid black line). The upper plot is for = 0.04 
and the lower plot for = 0.16. For comparison we show the 
fitting formula from Eisenstein & Hu (dashed black line), 
and our modification to their formula (dotted black line). The 
horizontal blue line is the approximation —8/^ from Q], and 
the horizontal dashed red line is our numerical solution to 
equation ([6]). 



works much better for small neutrino masses than its 
predecessor used to, as well as for very massive neutrinos. 
For 0.15 < n,n < 0.8, Qb/Qm < 0.3,/^ < 0.3, z = 
and Ni, ~ 3 the accuracy of the fitting formula is 
quite high. Note that the formula works equally well 
for = 0. Our revised code can be downloaded from 
|http://zuserver2. star. ucl.ac.uk7~lahav/nu_matter .power. f 



1 Accuracy of CAMB is set to 0.3 %. 



We also explore the effect of the dark energy equa- 
tion of state (w) on the matter power spectrum. For 
a constant f^, and in the regime of interest 0.02 < k < 
0.15 h Mpc"^, w affects the growth factor i5 which affects 
the matter power spectrum by changing the amplitude, 
but this can be altered by taking the ratio where this 
effect is minimal, i.e. is independent of w (to 4% in the 
regime of interest). 

For a fixed value of the parameter combination ^l^h, 
the impact of massive neutrinos on the matter power 
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spectrum is controlled by the fractional contribution of 
neutrinos to the total mass density in the Universe, ie. 
fiy = The scale where the suppression of power 
from neutrino free-streaming sets in is controlled by the 
comoving Hubble radius at the time when the neutri- 
nos became non-relativistic, corresponding to a comoving 
wavenumber 

kfs = 0.10a„ h^. (4) 

Note the obvious degeneracy with Q^^h. This param- 
eter, in models with negligible baryon density, sets the 
scale of the Hubble radius at matter-radiation equality, 
and so it determines the scale at which the matter power 
spectrum bends over in the case of massless neutrinos. 
For realistic baryon densities, there is an additional de- 
pendence on f2b- There is a statement, sometimes found 
in the literature, that the power spectrum depends only 
on fi, and not the indepenent values of the neutrino and 
matter densities. This is only true in the case where 
rim ft. and rib are held fixed as varies, and only for 
k > 0.8 h Mpc~^, since the neutrino free-streaming also 
enters this picture. However, this can be altered by nor- 
malising k by kfs (free-streaming scale). As shown in 
Figure ([2|), is a function of the ratio and 
Moreover, we see that (AP/ P)/ tends to -8 only for 
small neutrino masses whereas it goes to -4.5 for — 0.2. 
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FIG. 2: The dependence of AP/P at z = on fi™ and 
fl^, illustrating, via the scaling by kfs, equation Q, that 
just the ratio f„ = is insufficient to fully parametrize 
AP/P. From top to bottom: /^=0.2 (dashed line), 2 models 
with /^=0.1(S1^=0.03, Sl™=0.3 (dash-dotted hne), ^2^=0. 02, 
r2m=0.2 (solid line)), /,^=0.03 (dotted line). We see that the 
ratio ^p- 1 tends to a constant only for k > 50fc/s, but not 
to a universal constant. 

Exploring the E&H approximation to CAMB in Fig- 
ure ([2]), we observe that their formula provides best re- 
sults for only 1 massive neutrino and 2 massless neutrinos 
with large /y, whereas there is considerably less power on 
small scales for 3 massive neutrinos. However with our 
modified formula, this effect is altered and the approxi- 
mation is now valid for 3 massive neutrinos of any mass. 
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FIG. 3: AP/P at 2: = for 1 massive, 2 massless neutrinos 
(top panel) and 3 massive neutrinos (bottom panel), using 
CAMB (fuU line), E&H (dotted line), and our modified fitting 
formula (dashed line) .In all cases — 0.04, fim ~ 0.25, Qt = 
0.04. 

III. LINEAR GROWTH: AN ANALYTICAL 
APPROXIMATION 

The equation for linear evolution of density perturba- 
tions is 

6 + 2-S = 4TrGpa5, (5) 
a 

where 5 = Spm/pm, and pm and Spm is the density and 
the overdensity of matter, respectively. Light, massive 
neutrinos inhibit structure formation on small scales be- 
cause they free-stream out of the dark matter potential 
wells. Roughly, this can be taken into account by by mul- 
tiplying the driving term on the right-hand side of equa- 
tion ® by a factor ^ = i _ ^ n^. 

This equation is valid for k > 0.2 h Mpc^^. 

6 + 2-6^ AnGpoil- U)S. (6) 
a 

A. Einstein De Sitter Universe 

For an Einstein de Sitter universe,the solution to equa- 
tion (HD is given by 0, [13, [l| 
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S oc a^, 



(7) 



where p w 1 — Some authors (e.g. 0), have used 
this relation to estimate crudely the suppression of the 
power spectrum due to massive neutrinos ^ to be « 
— 8/,y. However, as we showed in the previous section 
this is a poor approximation. Our Figure [1] shows the 
limitation of this approach as it is only valid for very 
large scales fc > 0.5 hMpc~^ and small neutrino masses 
fi, < 0.05. More importantly, we have shown in Figure 
[5] that the suppression is not just a function of the ratio 
/j/, but of the matter and neutrino densities separately. 



B. A Dominated Universe 

We now evaluate the linear growth of perturbation, 
equation ([6]), for a Universe with a cosmological constant. 
On length scales below the present horizon dark energy 
does not cluster, and so it only affects H = a/a. For a flat 
cosmology with a dark energy component with constant 
equation of state p = wp, the Hubble parameter as a 
function of redshift z is given by 
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= - = Jfo Va„o(l + Z)3 + (1 - 0„,o)(l + z)3(l+»), 

a * 

(10) 

where Hq = lOO/i kms ^ Mpc is the present value of 
the Hubble parameter, parametrized by the dimension- 
less Hubble parameter h, and fimo is the present value 
of the matter density in units of the critical density that 
gives a spatially flat universe. The linear growth factor 
/ is defined by 



din (5 
din a' 



(11) 



^ The qualitative derivation for the suppression of P{k) goes as 
follows. From equation (0, the growth from matter-radiation 
equality epoch aeq to the present ao is 

S{ao) 



S{aeq) 



(1 + Zeq){l + Zeq)' 



(1- 



|f„ In(l-I-Zcq) 



(8) 



The power spectrum P(k) is the variance of the fluctutions <5 
in Fourier space, so massive neutrinos suppress it by the same 
factor as it suppresses S^, i.e. 



Pik,U)-Pik,U=0) 

p{k, u = 0) 
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fv ln(l + 2cq). 



(9) 



Conceptually this derivation contrasts two scenarios (with and 
without massive neutrinos) which yield at the present epoch the 
same amplitude of fluctuations. It also assumes that is the 
same for both scenarios, hence (1 + z^q) 23900n„h2 Yov 
the concordance model Q,mh^ = 0.175. This gives for the RHS 
of equation —9.6/^. Different coefficients may be 

obtained by taking into account whether the neutrinos became 
non-relativistic before or after matter-radiation equality. 
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FIG. 4: Density fluctuation 5{z) on comoving 8/i~^Mpc scale, 
normalized to CMB (top panel) and normalized to the value 
at 2; = (middle panel). The models shown in the flgure 
have {w,fv) equal to (—1,0) (solid line), ( — 1,0.04) (dotted 
line), (-0.8,0) (dashed line), and (-0.8,0.04) (dash-dotted 
line). In the case oi w = —0.8, we included dark energy 
perturbations according to CMBFAST. We have assumed a 
spatially flat universe, adiabatic fluctuations, and flxed the 
matter density 57m = 0.25, baryon density Q.f, — 0.04, the 
Hubble constant h = 0.7, scalar spectral index Us = 1, and 
the optical depth to reionization r = 0. The lower plot shows 
the density fluctuation S{z) as a function of at z = (solid 
line), z = 1 (dotted line), and z = 2(dashed line). 



and equation ([6]) can be rewritten in terms of / as 



df 
din a 



(12) 
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where a = (1 



is the scale factor of the universe, 



and fli-n{z) — H^fl^o{l + z)^/H^{z) is the time depen- 
dent density parameter of matter. The Runge-Kutta in- 
tegration of the set of equations pT|) . simultaneously 
gives the growth factor and the logarithmic derivative of 
the growth factor. Equations (fTT|) . (fT2|) are valid even 
when w evolves with time. 

Equation ([S]) is a simplified description of the effect 
of massive neutrinos on the growth of structures, but 
its solution is in good agreement with the results of 
detailed calculations with e.g. CAMB [l^. We have 
checked this by comparing In 5 found by solving ([T2|) to 
the value obtained from CAMB (expressed there as ag, 
the root-mean-square mass fluctuation in spheres of ra- 
dius 8 h~^Mpc). We note that the difference between 
the exact result and the solution of equation is at 
most 7.5 %. 

Since the simple prescription employed above to de- 
scribe the effect of neutrino masses on linear growth 
seems to work well over a range of parameters and cosmic 
epochs, we can motivate a simple analytical approxima- 
tion to the linear growth factor /. In matter-dominated 
cosmolo gies this quantity is well approximated by / « 
^m^{z) [ITI - fisj . An analytical approximation is also 
available for models with a cosmological constant, see 
e.g. p^ . For models without massive neutrinos but with 
a more_general dark energy component, Wang and Stein- 
hardt ^20] found that / = fi,'^,(-z), with 



where 



and 



ao 



5- 



w ' 

1—w 



ai 



3 (l-w)(l-3w/2) 
125 (1 -6u;/5)3 



(13) 



(14) 



(15) 



for a constant w. We should note that a is a very weak 
function of redshift for z > 1. 

To include the effect of massive neutrinos in the ap- 
proximation for /, we note that for massive neutrinos 
the growth is not only redshift dependent but scale de- 
pendent as well. By using CAMB, we calculated P(fc, z) 
over a range of values for f^, and varying w and ilm in 
the proximity of w = — 1 and — 0.25, we find that a 
reasonable fit for a flat universG, ^jyi-^^jj^ — 1, is given 

by 



(16) 



where 



fiik,u,nDE) - i-A{k)nDElu+B{k)fl~c{k)fl (17) 

Numerical values for A, B,C as a function of k are 
given in Table II. They vary over the redshift range z < 



k (scales) 


A(k) 


B(k) 


C(k) 


/ 


0.001 











0.825265 


0.01 


0.132 


1.62 


7.13 


0.824272 


0.05 


0.613 


5.59 


21.13 


0.815508 


0.07 


0.733 


6.0 


21.45 


0.811596 


0.1 


0.786 


5.09 


15.5 


0.806680 


0.5 


0.813 


0.803 


-0.844 


0.789606 



TABLE II: Numerical values for the coefficients A , B, C as 
a function of k from equation (|17|l . f for this values is given 
in the last column, for = 0.08, which corresponds to ap- 
proximately Sm^ = leV. For k > 0.5^ Mpc~^, fi is better 
approximated using equation p9p . 



10 up to 1%. Specifically for 8 h^^ Mpc we recorded the 
as given by CMBFAST, and the resulting approximation 
is given as: 

f{z;U,w,nDE) ~ {i~omDEU+i-9f^-9.8f^)nz,{z). 

(18) 

This formula is valid for < 0.15, and its accuracy is 
quite high for w = —l,z— 1, see Figure [5] Performance 
at w = —0.5 or at z = 10 is at most 2% worse than the 
w = — 1, z = 1 case. 

As expected /i = 1 on very large scales where k < knr- 
On very small scales, where k ^ knr, we would expect 
that / « [ri„i(l — /i/)]"° based on equation © and ([E 
ie. 



M = (1 - UY 



(19) 



Another approximation on these small scales follows 
from equation and by 2l|. However we note that 
this is based on the solution equation (O, which is only 
valid for an Einstein-de Sitter universe. We find that by 
solving the set of equations (fTT|) , for a A dominated 
universe, fi — I — 0.619/^, for z = 1, = 0.25, w — — 1, 
compared to ^ = p = 1 — O.6OO/1, of equations (l2]),([7l). We 
emphasize again this is only valid on very small scales, 
and for intermediate scales one should use equation (|17p . 



IV. THE w-f„ DEGENERACY 

Hannestad [l^] pointed out that cosmological neutrino 
mass limits are considerably weakend if one allows for 
w < — 1 (a case which is peculiar, as in this case the den- 
sity will increase with the expansion of the universe) . In 
his analysis, the parameter combination fl,^oh^ instead of 
fi, was varied, and the explanation of the degeneracy he 
gave was as follows: since f,y — rj^o/^^mo determines the 
small-scale suppression of the matter power spectrum, 
one can compensate for a larger Qi,q by increasing r^nio- 
If one assumes a fixed w — — 1 in the analysis, the Hubble 
diagram from supernovae Type la rule out values of f2niO 
much larger than 0.3. However, if one allows for w < —1, 
then the supernova data are compatible with consider- 
ably larger values of rimOi a-nd hence a higher value of 
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f -eqr (18) 

w=-0.8,z=1 
w=-1.0,z=1 
w=-0.8,z=2 
w=-1.0,z=2 



FIG. 5: Variation of the logarithmic growth factor / — ^ |" * 
with = and with redshift at a fixed scale 8/i~^Mpc. 
The ratio is shown as a function of for z=l and z—2. 

he following examples are shown: w = — 0.8, z — l(dashed 
line), w — -~l,z — l(dotted line), w — — 0.8, z — 2(dash- 
dotted line) and w — —1,2 — 2(dash-triple dotted line). The 
solid line is (1 — fu)°'°, from equation (|19p . and the long- 
dashed line is the full polynomial found in equation (|18|l , eval- 
uated at z = 1. Other parameters kept fixed Q.m = 0.3, Qb = 
0.04, /i = 0.7. 




FIG. 6: Contour plot showing the degeneracy of and Qde 
for ui = — 1 (full line) and w — —0.8 (dotted line), h — 0.7, 
Qb = 0.04, at z = 1 on 8 h~^Mpc scale. The contour levels 
correspond to values of / using equation (fT8)l . 



f2^o can be accomodated by a given value of f^. The de- 
generacy described above between fii^o f^nd w is indirect. 
We also emphasize again and again that the suppression 
is not just a function of f^,. 

The degeneracy is most relevant when the shape of 
the matter power spectrum P{k) at z = (or another 
fixed redshift) is used to constrain the neutrino mass 
(i.e. the bias between the galaxy distribution and total 
mass distribution is assumed constant and marginahzed 
over). Dark energy with a constant equation of state 
does not cluster on the scales probed by galaxy redshift 
surveys and does not affect the shape of the matter 



power spectrum, and hence w affects the limit indirectly 
through the mechanism described above. Furthermore, 
the degeneracy is almost completely broken when using 
the WMAP 3-year data as opposed to the first year Q. 

Based on the considerations in the previous subsections 
of this paper we would expect there to be a further and 
more direct degeneracy between and w: the same lin- 
ear growth rate can result from a range of combinations 
of values of w and In physical terms, this degeneracy 
can be understood as follows: neutrino free-streaming 
will supress growth of structure on small scales. How- 
ever, decreasing w (for fixed flmo) prolongs the era of 
structure formation, and reduces the value of the Hubble 
parameter in the matter-dominated phase. The Hubble 
parameter acts as a 'friction term' in equation ([5]), so re- 
ducing it will enhance the linear growth factor, see equa- 
tion (jlOp . To sum up, one can at least partly compensate 
for increasing by decreasing w. This degeneracy would 
be relevant if one were to use the growth of structure 
to constrain f^, and w. In Figure Q we illustrate this 
by plotting the root-mean-square mass fluctuation am- 
plitude S{z) as a function of redshift, for four different 
combinations of w and f,y. The middle panel shows the 
situation when if only the growth rate is measured (nor- 
malized at z = 0): in that case distinguishing between 
the four cases will require very accurate measurements. 
The situation is better when the absolute values of the 
fluctuations are measured, as shown in the top panel of 
Figure ([4]). This can be understood by noting that the 
absolute values depend, in the case when we normalize 
to the CMB on large scales, on w, but is only weakly 
dependent on f,y. The reason for this is that dark energy 
fluctuations are relevant on scales of the size of the hori- 
zon, and dominate integrated Sachs- Wolfe effect on the 
very largest scales. In contrast, the main effect of neutri- 
nos on the CMB is on smaller scales through a small shift 
in the position of the peaks and a slight enhancement of 
their amplitude. Thus, if the large-scale normalization 
is combined with a measurement of the growth rate, the 
degeneracy between w and fi, may to a large extent be 
broken. Figure [S] also shows the degeneracy of and w 
through the linear growth factor /. 

As indicated in Table [H there is an interplay in the 
roles of neutrino mass and the DE equation of state, 
when estimated observationally, which depend on geom- 
etry, the growth rate of perturbations and the shape of 
the power spectrum. We discuss briefly two probes where 
results from this paper, in particular the growth rate / 
in equation , could help in understanding degeneracy. 

The first is peculiar velocities (see e.g. [1^ for review). 
The rms bulk fiow is predicted in linear theory as: 

{v\R^)) = (27r2)-i ij2 J dkfP{k)W^{kR,) (20) 

where H/gC^-R*) is a window function, e.g. W{kR^) = 
exTp{—k^R^/2) for a Gaussian sphere of radius i?*. The 
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FIG. 7: The product D'^{f-l)P in equation (pTjl . using equa- 
tion (|18|) on 8/i~^Mpc scale, at £ = 50 (lower set of curves) 
and at £ — 100 (upper set of curves) for spatially flat universe 
models with Sim = 0.3, Qh = 0.04, h — 0.7, = 1.0, and 
w = -1, fi„ = 0.0 (fuU lines), w = -1, n„ = 0.01 (dotted 
lines), w — —0.8, flu = (dashed lines). 

velocity field at low redshift is insensitive to geometry. 
The power spectrum P{k) depends on the neutrino mass, 
but not on dark energy. Massive neutrinos would sup- 
press bulk flows [l^. However in any dependence of 
/ on neutrino mass and dark energy was ignored. Our 
equation (jlSp shows such dependence, which would re- 
sult in extra suppression of the bulk flows and in some 
degeneracy with w. This is also relevant on redshift dis- 
tortion. 

A second cosmological probe which depends on / is 
the Integrated Sachs Wolfe effect derived from cross- 
correlation of the CMB with galaxy samples [l^ . In the 
small angle approximation the predicted spherical har- 
monic amplitudes are (e.g. [26l[27j): 

(21) 



where Q{r) is a radial selection function for the galaxy 
survey, bg is the galaxy biasing factor, and D{t) is the 
linear theory growth function. The approximate relations 
for /, D and P{k) are useful to explore the degeneracies 
between neutrino mass and dark energy. As a simple 
illustration, we show in Figure ^ the product D^{f — 
1)P in the integrand above. We see that the selection 
function must extend to higher redshifts to distinguish 
between neutrino masses and w —1. See (28| for further 
analysis on the ISW effect in the presence of massive 
neutrinos. 

V. CONCLUSIONS 

We have demonstrated that the popular heuristic for- 
mula for the linear theory suppression of the matter 
fluctuations by free-streaming neutrinos, AP{k) / P{k) « 
— 8/i/, is valid only on very small scales (fc > 0.8/i/Mpc). 
However, it is not of practical use as this is in the strongly 
non-linear regime of matter clustering. We also provide 
a modified code for the E&H fitting formula. The linear 
growth factor in models with both massive neutrinos and 
a dark energy component with equation of state param- 
eter w =constant, has been calculated numerically, and 
we have provided a fitting formula justifled by simple ar- 
guments. Furthermore, we have explained the indirect 
degeneracy between w and rriiy found in [l^ , and argued 
that there is a further, more direct degeneracy between 
w and fi, when the linear growth rate of density pertur- 
bations is considered. This degeneracy can, however, be 
lifted by measurements of the absolute scale of the mass 
fluctuations, for example from the CMB. 
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